NASA /TM— 1999-209082 



The Reconstruction Problem Revisited 


Ambady Suresh 

Dynacs Engineering Corporation, Inc., Brook Park, Ohio 


National Aeronautics and 
Space Administration 


Glenn Research Center 


May 1999 


Acknowledgments 


The author would like to thank H. T. Huynh and D. Paxson for several interesting discussions on upwind schemes. 
The author is also grateful to the Computing and Interdisciplinary Systems Office at 
Glenn Research Center for permission to publish this work. 


Trade names or manufacturers' names are used in this report for 
identification only. This usage does not constitute an official 
endorsement, either expressed or implied, by the National 
Aeronautics and Space Administration. 


Available from 


NASA Center for Aerospace Information 
7121 Standard Drive 
Hanover, MD 21076 
Price Code: A03 


National Technical Information Service 
5285 Port Royal Road 
Springfield, VA 22100 
Price Code: A03 
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Ambacly Suresli 
Dynacs Engineering, Inc., 
Cleveland, Ohio 44135, USA. 


Abstract 

The role of reconstruction in avoiding oscillations in upwind schemes is reexamined, 
with the aim of providing simple, concise proofs. In one dimension, it is shown that if the 
reconstruction is any arbitrary function bounded by neighboring cell aveiages and incieasing 
within a cell for increasing data, the resulting scheme is monotonicity preserving, even t hough 
the reconstructed function may have overshoots and undershoots at the cell edges and is in 
general not a monotone function. In the special case of linear reconstruction, it is shown 
that merely bounding the reconstruction between neighboring cell averages is sufficient to 
obtain a monotonicity preserving scheme. 

In two dimensions, it is shown that some ID TV I) limiters applied in each direction 
result in schemes that are not positivity preserving, i.e. do not give positive updates when 
the data are positive. A simple proof is given to show that if the reconstruction inside the cell 
is bounded by the neighboring cell averages (including corner neighbors), then the scheme 
is positivity preserving. A new limiter that enforces this condition but is not as dissipative 
as the Minmod limiter is also presented. 


Introduction 

In this paper, we begin by reexamining the reconstruction step in upwind schemes 1- *’ 
with the aim of deriving slightly more general results and providing simple, concise proofs. 
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To introduce the one dimensional problem, consider the model equation 


+ Ur — 0 , 


( 1 ) 


where / is time, x is distance, and = it 0 (*) is the initial condition. We denote by u 1 - 

the cell averages of u(xj n ) on a uniform grid of width A.r and cell center X; = j&x. 


1 f x 

A,/, 


j + At/2 
— Ax/2 


u(x, /" ) dx. 


( 2 ) 


Let R'j((,). £ = (x — Xj ) / A // denote the reconstruction of u(xj n ) in the cell j expressed 
in local coordinates. The new cell averages are obtained by converting the reconstructed 
profiles a time step and averaging the resulting profiles over a cell. The result is 

= lZ^- At)dt + !-m 

where A = A t/h. . For clarity, we omit the superscript n when there is no confusion, i.e <7 • 
denotes 

For increasing data, i.e. uj > uj-\ for all j . we are interested in schemes for which 
the new cell averages are also increasing. Such schemes are called monotonicity-preserving 
schemes. They mimic the exact solution and as a special case it follows that a step function 
propagates without spurious oscillations. 

The question we are interested in is: 

What art sufficient conditions on the reconstruction RjH) that will ensure that the scheme 
(3) is monotonicity-preserving Y 

This question has been studied extensively over the last two decades in the upwind liter- 
ature, and we review existing conditions below. Our aim is to provide a new, simple, concise 
proof that covers the general case. Generality is achieved here by considering a) completely 
arbitrary functions for reconstruction and b) by considering Monotonicity-Preserving (MP) 
schemes as opposed to Total Variation Diminishing (TVD) schemes. TVD schemes are 
monotonicity-preserving, but not vice-versa. In addition, MP schemes can be designed to 
avoid the chronic loss of accuracy at smooth extrema incurred by TVD schemes (see Kef. 7 
for an example). 


Existing Conditions 

There are two different conditions stated in the literature. The first condition, introduced 
by Van Leer 1 in his seminal paper on upwind schemes and often called Van Leer’s condition. 


NASA/TM— 1 999-209082 


2 




Figure 1: An oscillation created by reconstructions bounded by neighboring cell averages. 
A = 0.25. 

states that the reconstruction RjU) must be bounded by the neighboring cell averages, i.e. 
for increasing data., the condition is 

iij - 1 < m) < «i+! (1) 

While this condition suffices for linear reconstructions, we will show by an explicit coun- 
terexample that this is not sufficient for general reconstructions. Indeed, let «j _ 3 = -6. 
«._ 2 = -3. 0 J _ 1 = 0, u 3 - 1, Ttj+ 1 = 4, and the reconstructions be Rj- 2 (Z) = -3 + 6£; 

= (4 - 24£ - 4S£ 2 )/7; /?,(£) = 1 - (4 - 24/ - 48/ 2 )/7; R j+ 1(0 = 4. These re- 
constructions are plotted in Fig. 1. Note that the data is increasing and (4) is satisfied in 
each cell. However, the cell averages at the next time step given by (3) for A = 1/4 are 
tf}+l = 27/112, tt” +1 = -3/28, and = 103/28. So monotonicity is not preserved. 

Another condition (Harten et. al 2 , LeVeque 3 ) is stated in terms of the total variation of 
the reconstruction. Let TV{f(x)) denote the total variation of /(. r) given by 

TV(f)= / +OC \f'(x)\dr (5) 

J — co 

and TV {it) denote the discrete total variation of the cell averages given by 

TV(u)= £ I'L+i - "/I (6) 

/.■= — oo 

Then the required condition on the reconstruction is that 

TV(R{.)) < TV {it). (7) 
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This condition holds for any general reconstruction. For increasing data, (7) implies that R{.) 
is an increasing function with no overshoots at the cell faces. However, designing schemes 
that satisfy this condition is a bit more difficult because the conditions on Rff) now involve 
Rj- i(£) an( l Rj+ i(£) hi addition to the cell averages. This condition is also unnecessarily 
restrictive since most of the well known TVD schemes do not satisfy it. 

For example, consider the data given by = 0, iij = 0.2, « i+1 = 0.8, and u j+2 = 2.6. 
The reconstructions obtained by the popular Superbee limiter and Van Leer's Average limiter 
( (15) and (46) in the appendix) are identical and given by Rj(£) = 0.2 + 0. 4£ and R j+ j(£) = 
0.8 + 1.2 £. which has an overshoot at the face x J+1 / 2 . Thus both these schemes violate (7) 
although both limiters are monotonicity preserving. This example clearly illustrates that 
condition (7) is unnecessarily restrictive for monotonicity preservation. 

Thus, ol the two conditions available in the literature, neither are completely satisfactory. 
Van Leer's condition fails for general reconstructions while the second condition (7) is so 
restrictive that many well known schemes do not satisfy it. In the next section, we introduce 
a slight modification of Van Leer s condition that is more general than (7) and for which 
monotonicity preservation can be proved for general Rj{£). 

A Modified Condition 

First, we need two inequalities that hold for increasing functions. A function f(x) is an 
increasing Junction i( for any x ] > x 2 , j ( x i ) > / ( x 2 ) ■ Let f ( x ) be a .11 increasing function in 
[—1/2, 1/2] and let its mean over this interval be denoted by /. Then, for 0 < A < 1, we 
have the following two inequalities that are useful. 

1 f 1 / 2 - ] /-1/2— A 

A L-> m)d( - f ' —xL I2 tm</- (8) 

The first inequality states that sampling an increasing function from the right yields an 
average greater than the cell average. To prove this inequality, note that if we denote by 
fj( A) the left hand side of the first inequality, we get A</'(A) = /(1/2 — A) — </(A). As this is 
nonpositive for all A, the inequality follows. The proof of the second inequality is similar. 
The modified condition can be stated as follows: 

Theorem-1: Assume we have increasing data and the reconstructions Rj(£) satisfy the three 
conditions give n below for all j : 

mi =*> (9) 

»./-i < m) <«t+i (io) 

Rj{0 is an increasing f unction (11) 
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Then, for 0 < A < 1, the cell averages at the next time step are also increasing and (3) 
becomes a monotonicity preserving scheme. 

Proof: 

Since Rj- 1 (£) is an increasing function, the first integral in (3) satisfies 

{ / 1/2 Rj-iit)# > '0-1 ( 12 ) 

A J 1/2- A 

and from (10) the second integral satisfies 

/ ,/2 ' A R,U)d{ > (1-A )«j-i 03) 

Thus (7" +) > tij _ j . Similarly, from the bounds on f?j_ 2 (£)* we have 

/ 1/2 Rj- 2 (Odf < A (14) 

J 1/2— A 

and since F?j_i(£) is an increasing function, 

-p— — r I'' 2 ’" Rj-i(f)df < uj (15) 

1 — A .7-1/2 

So < Uj_i. Hence, f<" +1 > completing the proof. D 

Thus, if we add to Van Leer’s condition the requirement that the reconstruction be 
increasing when the data are increasing, we can prove monotonicity preservation for any 
genera] reconstruction. Note that conditions (9), (10) and (11) are more general than . (7) 
since TV{R{.)) can be greater than TV(u). In addition, all the well known second order 
TVD schemes satisfy (9), (10) and (11) so that Theorem -1 is not unusually restrictive. 

Although adding the requirement that the reconstruction be a monotone function is a 
fairly trivial modification, we have not seen Theorem - 1 stated and proved anywhere in the 
literature. Indeed, this is all the more surprising since the fairly restrictive condition (7) is 
widely quoted. Monotonicity of the reconstruction was imposed as a condition by Cololla 
and Woodward 8 in the desigii of their piecewise parabolic method (PPM), but monotonicity 
preservation is not proved there. 

TVD Schemes 

TVD schemes have some theoretical advantages over monotonicity preserving schemes 
such as a guaranteed convergence to a weak solution of the underlying conservation law 
T Monotonicity-preserving schemes differ from FY D schemes only in theii tioatment of 
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local extrema. It is not surprising, therefore, that with some recipe for reconstruction at an 
extrema, a result similar to Theorem - 1 can he proved for TVD schemes. 

Theorem-2: Let the reconstructions Rj(f) satisfy the conditions given below for all j : 




m) = » 

i! 

(16) 


tv hi 

m uj-\ < uj 

+ 

VI 


< m) 

• ^ 

VI 

and Rj(() it 

i an increasing f unction ; 

(17) 


iv hi 

:U Uj-l > Uj 

> Uj + l ■ 


Uj-\ > Rjtt) 

+ 

* 

Al 

and Rj(Z) i 

s a decreasing f unction ; 

(18) 



otherwise 

: 




3? 

cn 

II 

’./* 

(19) 


Then . for 0 < A < 1. (3) becomes a TVD scheme. 

Proof: We introduce some terminology on data types. Let us say the data is increasing at 
j when it ,- 1 < Uj < iij+\ and decreasing at j when Uj_i > uj > uj+i- From Theorem- 1 it 
lollows that for data increasing at j and j — 1 , or decreasing at j and j — 1 the new cell 
averages lie inside the range defined by neighboring old cell averages, i.e. 

Min(uj,«j_i) < i/" +1 < (20) 

Consider the case where data is increasing at j but we have a local minimum at j — 1, i.e. 
Uj — i < Uj-i and u,-\ < a , . In this case using (19) and (3) we can show that the bound (20) 
is still satisfied. Similarly, we can show that the bound (20) holds also for all other cases 
such as data decreasing at j — 1 and a minimum at j. or increasing at j — \ and a maximum 
at j etc. In short, (20) holds for all data. 

Since (20) holds for all data, no new extrema can be created and the value of a local 
maximum can only decrease and the value of a local minimum can only increase. Hence the 
total variation can only decrease and the theorem is proved. □ 

The Linear Case 

The linear case has been exhaustively studied ,-<i in the context of TVD flux limiters 
and will not be repeated here. However, we wish to highlight the following peculiarity 
of the linear case, namely that monotonicity is preserved even without requiring that the 
reconstruction be increasing for increasing data, i.e. even without (11). In other words, for 
increasing data, the slopes can be decreasing and monotonicity is still preserved so long as 
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the linear reconstructions remain bounded by neighboring cell averages. For completeness, 
a proof is sketched out below. 

Assume increasing data and linear reconstructions i?j(£) = Uj + satisfying (10) loi 
all j. The case where u } = w,_i is trivial. For Uj ± iij-i, by a suitable normalization we can 
take the data to be u J+ 1 > 1, fij = T Uj- \ — 0 arid i/j _ 2 £ 0* A direct calculation of (3) 
gives 

u" +1 — == 1 — A — Ai/j _2 + A(1 — A)(2.s y _i — «sj — Sj-'i)! £ (^1 ) 

From (10), |^-| < 2, |$j-i| < Min( — 2t7 ? _ 2 , 2), |^_ 2 | < -2«j- 2 . Under these bounds, it can 
be verified that the right hand side ol (21) is positive lor 0 < A < 1 and u ,-2 < 0. Thus 
monotonicity is preserved. 

Note however that without condition ( 11 ), ■«* +1 can lie outside of iij]. So here we 

have an instance of the scheme being monotonicity preserving without being bounded by 
the initial data. 


Reconstruction in Two Dimensions 

For the two-dimensional advection equation 

u t + au x + buy = 0 (22) 

with initial condition u{x,yA)) = u 0 (x,y), a. set of conditions under which higher order 
reconstructions preserve monotonicit.y is not known. 1 lie focus then shifts to schemes that 
are positivity-preserving (PP) 9-14 where some theoretical results can be proved. These are 
schemes that give positive updates when the initial data are positive, or equivalently, where 
the cell averages at time n + 1 are bounded by the cell averages at time level n. In the 
meteorology literature, such schemes are also referred to as positive definite. 

In one dimension, the boundedness result used to prove Theorem- 1, i.e. that 

Uj - 1 < uf' < uj ( 23 ) 

suffices to prove monotonicity. In two dimensions, this is no longer the case and PP schemes 
will not preserve monotonicity in general. 

To introduce the reconstruction problem here, let us begin with a base scheme that uses 
midpoint rule in time and linear reconstructions in each cell. This scheme is chosen since: 
a) it can be readily extended to the Euler equations, b) it requires only one Riemann solver 
calculation per lace per timestep and c) it has already been implemented in commeicial 
codes 15 . 
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We assume a uniform grid of spacing Ax and Ay and a > 0 and b > 0. A linear 
reconstruction is assumed in each cell that can be written in normalized coordinates as 


Ru = + Sfjt + sfjt, 


(24) 


with |£|,|r/| < 1 /2. The scheme can be written as 


u n + l 


= 


K(' 


n+l/2 

1+1/2J 


,. n + 1 / 2 v _ , / n+l/2 

U i-\/2,j) A y\ U i,j+ 1/2 


— U 


n+l/2 v 
ij — 1/2 J 


(25) 


where A, = a At /Ax, X y = bAt/Ax and the interface values are obtained by taking a half 
timestep inside the cell and can be derived as 


+ 1/2 
*+1/2, J 
+ 1/2 
* 0 + 1/2 


= U n ■ 
W «0 


*o 


+ (l-+)A',y2- KSfJ‘2 

+ (1 - A ,)Sy-2 - X.Sfj/2 


(26) 


The linear scheme with Sfj and Sfj given by their central difference values such as 


Sfj = (u ( - +w -« i . 1 j)/2 (27) 

is second order accurate in time and space and stable in the region A* > 0, A a > 0, A^+Ay < 1. 
Thus the time step can be defined in terms of the CFL number by 


At = CFL/ (a/ Ax + bj Ay) (28) 

We can now state the reconstruction problem in two dimensions: 

What are sufficient conditions on the reconstruct ion //,,(£, r/) that will ensure that the 
scheme (25) is positivity preserving ? 

For example, if we use any of the TVD limiters to calculate the x and y slopes, is the 
resulting scheme positivity preserving ? The answer to this depends on the limiter as the 
following example shows. Consider negative data given by 


< 1 - 1,1 — — 1 • t<o,i — 0 

11 - 2.0 = — 1 00 «_i,o = — 1, iio.o = 0 a i,o = 0 

M-1,-1 = -1 «o,-i = -10, Ui,_ x = -100 

11 o ,-2 = — 100 

For Aj. = 0.6, A v = 0.2, both the Superbee (46) and the Average limiter (45) (these 
limiters are defined in the appendix below) give = 16/50, which shows that these 

schemes do not preserve positivity. The Minmod limiter, however, gives Uqq 1 = —57/50, 
which is acceptable. A closer look at the reconstructions in this example reveals that for the 
Superbee and Average limiters, the reconstructions in each cell fall outside the range of the 
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cell averages of its neighbors. This condition turns out to be sufficient to prove positivity as 
shown below. 

Let us define Nij to be the set of cell averages of the immediate first, order neighbors ol 
the cell (i.j) and u tJ . 


At, = 




(30) 


and Uij to be the range of variation of u on A (J . 


‘..i 


[Min [Ay,,], Max [Ay,.,]] 


(31) 


Theorem-3: Assume that for the scheme (25) the reconstruction in each cell is bounded by 
its immediate neighbors, i.e. for all (i,j) 


H 


ej 


e 1 


ij 


(32) 


Then , for X x + A v < 1 


-n+1 


€ t ij U L 'i- \J U [ i,j- 1 


i.c. the cell average at the next time step lies inside the union of the averages in the neigh 
borhoods of (> — Lj) and (i , j — 1). 

Proof: Using (25) and (26) , we can write the scheme as 

= (1 - A, - A y )(f/i,j - K.%/2 - Xy Slj/2) 

+\ r (ui-i,j + (1 A r ) <h;_ip/2 — A y 2) (34) 

+Xy(u i , J - l ~ X f Sfj-J 2 + ( 1 — A y ) S?j_ i / 2) 


which can be written in terms of the normalized reconstructions as 

= (1 - A,. - A v ) Ri.j{—Xj:/ 2, —X u / 2) 

+X x - X T )/ 2,-Xy/ 2) (35) 

+A v f?,,_ 1 (-A,/2.(1 - A„)/2) 


Under the stability limit A r + X y < 1, the points where the reconstructions are sam|)led in 
the above equation lie inside their respective cells. Thus the right hand side is a convex 
combination of sampled reconstruction values which are bounded by (32) and thus (33) is 

immediate. ^ 

It is easily verified that Theorem-3 holds true for all propagation directions provided 
the interlace values (26) are defined from upwind cells. 1 he geneial ( FL limit is then 
|Ar| T | Ay | < I- 
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The condition that the reconstruction be bounded over the whole cell by all its immediate 
neighbors (including corner neighbors) was introduced by Barth and Jesperson 16 in the 
context of unstructured grids. The above result shows that this condition is useful for 
constructing PP schemes in multi-dimensions for structured meshes as well. 

We remark that the choice of stencil, i.e. (30) used above is somewhat arbitrary. In fact, 
for the Minmod scheme, the reconstruction is so tightly constrained that R t j actually lies 
inside the range of tfij, and two linearly independent cells from tq± jj, and u,j±i . From this 
it follows that is bounded by the four values u. hJ , ii t _ X j, u l<3 -\ , and , a much 

smaller stencil than (33). 

In the next section, we derive an efficient limiter that enforces (32) but is not as dissipative 
as the Minmod limiter. 


A PP Limiter 


The procedure for modifying the slopes to satisfy (32) is far from unique. The approach 
described here is similar to the approach of Barth and Jesperson 16 , adapted to structured 
grids and modified so as to depend continuously on the data. 

To satisfy (32) within a cell, the idea is to restrict the slopes so that the reconstructed 
values at the four corners of the cell lie inside the required interval, i.e., (31). Due to the 
symmetry of the grid, this is equivalent to restricting |5^| + | S v t J | to lie inside another 
interval. We skip the details and give the result in algorithm form. 

Let the slopes be defined initially by their centered difference values, namely 


5'? . = 

s y = 

hj 


t/i j_i)/2 


(36) 


and let \ min and V max be defined bv 


1 min — Mm 


mar — MflX 


+ l UiJ, ^2 + 1, i+1 Utyj' 

Mi — lji 1 ,j 

^f — 1J— 1 'Mjjii t-Hj — 1 *b + lj — 1 i 

n*—! j-j-i i j 4- 1 ^ v i r j -e 1 ^i,j ? 


(37) 


Ui-lJ - Uij , 


e, 


u 


*+ I,i 


i i,j 


Ui — lJ — l UiJ'f ^ij— 1 + — 1 U t j 


where t is a small positive number (< = 10 10 in all our numerical experiments). Then we 
define 

v = 2 Min(|V^„|,|V^ ...h 


The final slopes are then 


\Sfj\ + |S?J) 

$■, = Min(l , V) Sfj 
Sfj = Min(l , V) Sfj 


(39) 


(40) 
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It can be verified that after limiting, the reconstructed values in the cell (Lj) lie inside the 
interval and that the reconstruction depends continuously on the data. 

The main problem with this and other limiters is the loss of accuracy at extrema which 
make these schemes only first order accurate in the Max. norm. If a scheme is required to be 
second order accurate everywhere (i.e. in the Max. norm), such a scheme will advect. second 
order polynomials exactly. It follows that while advecting smooth extrema such a scheme 
will give updates not bounded by the data. It thus appears that, second order accuracy in 
the Max. norm and positivity preservation are mutually exclusive. 

This is not to say that some compromise might be found w'here positivity preservation 
mav be relaxed in some useful way near extrema. However, such a compromise has eluded 
us thus far. 


Numerical Experiments 


We present some numerical experiments on two dimensional advection to assess the per- 
formance of PP schemes. We solve (22) on the domain [-1,1] X [-1.1]. Periodic boundary 
conditions are imposed at. the boundaries. A number ol initial conditions, advection diiec- 
tions and timcsteps were explored but only two cases are reported here. After computation, 
we scan the computed solution lor overshoots and undershoots and present line plots thiough 


these regions. These represent the worst case results. 

Example 1: In this example, a uniform grid ol 120 X 120 cells is used with convection 


velocities a — 0.8 and 
final time is t = 20, 
condition is 


h — 0.1. The timestep was calculated by (28) with CFL — 0.8 and the 
which corresponds to 8 periods in x and one period in y. llie initial 

u(x,y< 0) = 1 for ( x^ T ^ )* / 2 <0.4 (41) 

— 0 otherwise 


In Figure 2, the results along the line y = -0.264 (j=45) are showrn. The solid line is the 
exact solution. As can be seen both the Superbee and the Average limiter give significant 
overshoots while the Minmod and PP limiters have no overshoots or undershoots. The PP 
limiter is not as dissipative as the Minmod limiter but shock resolution is still quite poor. 
In terms of efficiency, the PP limiter is m the same ball park as any ot the T\T) schemes. 

Example-2: The initial condition here is 


u{x,y, 0) = exp( — + y 2 )) 


( 12 ) 


with 3 = 200. The convection velocities are a = 1 and 6=1. The same grid and CFL 
number are used as in Example- 1, but the final time is / = 2 which corresponds to oik- 
period in both .r and ij directions. The results along the line y — 0 are shown in higure b 
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where the loss of accuracy at extrema is clearly visible. The new' PP limiter gives results 
similar to the Average limiter. The results from the Minmod limiter are much poorer and 
are not show'll. 

We remark that a large number of numerical experiments with various initial conditions 
seems to suggest that all the TVD limiters are positivity preserving when the advection 
direction is 45 degrees. 

Conclusions 

In this paper, we have considered anew what conditions a reconstruction must satisfy for an 
upwind scheme to have desirable nonoscillatory properties. In one dimension, a simple con- 
cise proof is presented that if any general reconstruction is increasing for increasing data and 
bounded by neighboring cell averages, the resulting scheme is monotonicity-preserving. In 
the special case of linear reconstruction it is showm that merely bounding the reconstruction 
between neighboring cell averages is sufficient for monotonicity preservation. 

In two dimensions it is shown that if the reconstruction over the whole cell is bounded 
by the averages of all its first order neighbors, then the scheme is positivity preserving. An 
efficient limiter that, achieves this is also presented. 

We hope that these concise proofs will be useful in the teaching of upwind methods and 
also lead to the design of new' schemes. 


Appendix 


For completeness, w r e summarize several popular one dimensional TVD limiters here. The 
reconstruction in each cell is assumed to be Rj(£) — uj + Sj£ and these limiters are various 
recipes to determine nj from the data. Let s + and .s_ denote the forward and backward 
differences given by 

•S+ = (Uj + l - Uj) 

- s - ={Uj-Uj- 1 ) 

The Minmod limiter is given by 


- s j = ^ [sgn(5+) + sgn(s_ )] Min(|a_|, |.s + |) 

The Average limiter of Van Leer can be written as 

8 j - [ s S n (*+) + sgn(s-)] Min((s+ + s_)/ 2. 2|s_|, 2|s+|) 

and the Superbee limiter of Roe can be written as 


(43) 


(44) 


(45) 


- [sgn(s + ) + sgn(s_)] Min(Max(|.s + |, |s_|), 2|s_|, 2|.s H 


(46) 
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As mentioned above, the Minmod limiter is the most diffusive and the Superbee the least 
diffusive. 
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